Precise Control of Molecular Dynamics with a Femtosecond Frequency Comb - A 
Weak Field Route to Strong Field Coherent Control 
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We present a general and highly efficient scheme for performing narrow-band Raman transitions 
between molecular vibrational levels using a coherent train of weak pump-dump pairs of shaped 
ultrashort pulses. The use of weak pulses permits an analytic description within the framework of 
coherent control in the perturbative regime, while coherent accumulation of many pulse pairs enables 
near unity transfer efficiency with a high spectral selectivity, thus forming a powerful combination 
of pump-dump control schemes and the precision of the frequency comb. The concept is presented 
analytically and its feasibility and robustness are verified by simulations of dynamics in Rb2. We 
consider application of this concept to the formation of stable, deeply bound, ultracold molecules. 

PACS numbers: 33.80.-b, 33.80.Ps, 34.30.+h, 42.50.Hz, 42.65.Dr, 33.90.+h 



Although mode-locked lasers emit broadband ultra- 
short pulses, they can be utilized to perform frequency 
selective excitation just like narrow-band CW lasers due 
to their precise frequency comb 0. This spectral se- 
lectivity is explained by the very long inter-pulse phase 
coherence, which allows for coherent accumulation of the 
excitation amplitudes from multiple pulses in an excited 
material system, similar to a generalized Ramsey experi- 
ment. This idea led to the realization of direct frequency 
comb spectroscopy in atomic systems Q- Here we pro- 
pose to apply the principle of coherent accumulation, 
combined with weak field coherent control, to precisely 
control molecular dynamics at high efficiencies. 

While analysis of coherent quantum control is rela- 
tively simple in the weak field perturbative domain 
extension to strong fields is not straightforward. Analytic 
models exist only for simple cases 0,0 and solutions are 
often found by numerical optimizations 0. The core of 
our approach is to exploit analytic perturbative models 
to design "ideal" weak pulses and to achieve the high 
overall efficiency by coherently accumulating many such 
pulses. This avoids the complication of strong held de- 
sign while gaining high spectral selectivity offered by the 
frequency comb. Since the maximal number of accumu- 
lated pulses is inherently limited by the coherence time of 
the material ensemble, we expect our approach to be ap- 
plicable particularly well to ultracold atomic / molecular 
ensembles, where coherence times are long. 

Motivated by the goal to produce deeply bound ul- 
tracold polar molecules from loosely bound Feshbach 
molecules, we consider the Raman transition shown in 
Fig. ^ from a single vibrational level embedded in a 
dense environment of other levels near the dissociation 
limit, to a single deeply bound vibrational level. Loosely 
bound Feshbach molecules can be generated with high ef- 
ficiency via adiabatic sweeping through a magnetic Fesh- 
bach scattering resonance in an ultracold atomic ensem- 
ble Ml. Since molecules are not amenable to standard 



laser cooling techniques, magneto-photo-association of 
cold atoms is a major avenue pursued for obtaining cold 
molecules, which now represents one of the most exciting 
research fields in cold matter. So far, stable, ultracold 
ground state molecules have not been produced in high 
densities mainly because of the poor wavefunction 
overlap between the spatially spread continuum states of 
colliding atoms and the localized molecular states. Fesh- 
bach molecules appear therefore to be an important mid- 
stage towards deeply bound ultracold molecules. 

When the energy level spacing near the input state 
is small compared to the natural line-width of excited 
states, the application of CW techniques is not straight- 
forward. Any scheme relying on populating an excited 
state during the transition is inherently limited because 
most of the population will be lost to spontaneous decay 
during the time needed to resolve the input state. Short 
pump-dump pulses could eliminate spontaneous emission 
losses but then the total interaction time is too short. 
A similar problem is met by stimulated Raman adiabatic 
passage (STIRAP). For the adiabatic passage to hold, 
the Rabi frequency f2 must be larger than the decay 
rate T from the intermediate state: fi > T [|. Since 
both the power broadened pump line- width and the two- 
photon line-width of STIRAP are of order fl, the possi- 
bility of population transfer to neighboring levels should 
be considered. Thus, predictions for using STIRAP to 
photo-associate an atomic Bose-Einstein condensate de- 
pend strongly on modelingof the decay process and the 
atomic continuum states [lTJ . 

To overcome these difficulties, we employ the scheme 
illustrated in Fig. Q which is based on a phase-coherent 
train of shaped pump-dump pulse pairs. Each pulse pair 
is weak, i.e., it transfers only a small fraction of the in- 
put population to the target state. Coherent accumu- 
lation then enables a high overall transfer efficiency. In 
the time domain, each pump pulse excites a wave packet 
that starts to oscillate in the excited electronic potential. 
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FIG. 1: Basic Raman control scheme, (a) Typical electronic 
potentials and vibrational levels (here, Morse potential fits of 
RD2). Population is transferred from the input vibrational 
level near the dissociation limit (v 1 = 130) to a deeply bound 
level (v = 45), mediated via a broadband wavepacket in the 
excited electronic potential, (b) Time domain picture. A 
train of phase coherent pump-dump pulse pairs interacts with 
the molecule. Pulse pairs are shaped to achieve efficient pop- 
ulation transfer. The intra-pair time is half the vibration time 
(i~ V ib) of the intermediate wavepacket and the inter-pair time 
is the repetition time of the source l// rep . (c) Frequency do- 
main picture. A tooth-to-tooth match between the pump and 
the dump frequency combs, with f rep equal to a sub-harmonic 
of the net Raman energy difference, locks the relative phase 
between pulse pairs to the free evolving Raman phase. 



After half a vibration, this wave-packet reaches the in- 
ner turning point where the dump pulse drives it to the 
target state. Since population appears on the excited po- 
tential only for half a vibration, this scheme eliminates 
spontaneous emission losses. In between pulse-pairs the 
system is left in a coherent superposition of the input and 
target states that evolves freely. Therefore, in order to 
enable coherent accumulation of population at the target 
state, the temporal phase difference between pulse pairs 
Aip should match the phase of the free evolving Raman 



coherence. In the frequency domain, the combs of the 
pump and the dump pulses must overlap tooth to tooth 
and the repetition rate f rep must match a sub-harmonic 
of the Raman energy difference. 

Let us clarify the role of the different amplitudes and 
phases involved. All pulse pairs share a common tempo- 
ral (spectral) shape and phase (E(t) or E(lu)), which is 
designed to maximize the transfer efficiency for a single 
pulse pair. The relative phase between successive pulse 
pairs, Aip, is then controlled via stabilization of the fre- 
quency comb to achieve coherent accumulation. 

For the very first dump pulse to drive all the excited 
population to the empty target state its "pulse area" 
should be ~7r. After the second pump pulse however, 
the excited population is about equal to the population 
already in the target state, so now only a ir/2 "area" 
dump pulse of the appropriate phase is required to per- 
form the transfer, just like in a Ramsey experiment. In 
general, the fraction p of population excited (or dumped) 
is related to the "pulse area" A by sin 2 (A/2) = p. As 
population is accumulated in the target state and de- 
pleted from the input state, the dump "area" of the nth 
pulse Ad[n] should decrease and the pump "area" A p [n] 
increase for each pulse pair to transfer the same fraction 
of population (for N pulses), 
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The "pulse area" is not a well defined quantity outside 
the context of a two-level system. However, the ratio 
of populations in the excited wave-packet and the in- 
put (target) state defines an effective area for the pump 
(dump) pulses. This concept proves useful mainly for 
weak pulses, where the excitation is predominantly a one- 
photon process. In general, for any given pump series, a 
dump series can be matched according to the nth fraction 
of excited population p[n]. Clearly the very first dump 
pulses and the last pump pulses are of areas near tt and 
cannot be considered weak, but for a large N, the major- 
ity of the population is transferred by the accumulative 
effect of all pulses, which are mostly weak. 

For an efficient pump-dump process it is required that 
the wave-packet \ip p ), excited by the pump from the input 
state \i) and propagated for half a vibration, will overlap 
perfectly with the wave packet \tpj) that would have been 
excited from the target state \t), by the time reversed 
dump. For weak pulses we can express these two wave- 
packets, using first order perturbation theory, as 

\ip p ) = ^2E p (uj)exp[i(j) D (Lj)} \uj) (uj\d el \i) 

UJ 

|^>=EWM<*|*>, (2) 

UJ 

where \ui) denotes the vibrational states in the excited 
potential using the detuning to from the pulse carrier 
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FIG. 2: Franck-Condon factors ({v\v')) vs. frequency detun- 
ing from the central state of the wavepacket; (a) for excitation 
from level v' = 130 in the ground electronic potential to levels 
centered around v = 162 in the excited potential and (b) for 
de-excitation of these levels to v = 45 in the ground potential. 



frequency as a vibrational index, F p {uS) = (u)\d e i\i) and 
Fd(u>) — (u>\dei\t) are the pump and the dump transition 
dipole matrix elements, which under the Condon approx- 
imation are propotional to the Franck-Condon factors 
(uj\i) , (uj\t) . d e i is the electronic transition dipole mo- 
ment and Ep(E^) is the spectral amplitude of the pump 
(time-reversed dump) field. </>d(w) is the spectral phase 
acquired by the wave-packet between the pulses, which 
reflects both the delay of half a vibration and the disper- 
sion of the wave packet as it oscillates in the anharmonic 
excited potential. For the Morse potential fit, used later 
in the simulations, analytic expressions exist for the vi- 
brational states and energies yjj , so both the dipole ma- 
trix elements and </>_d(w) are known. 

As a result, perfect overlap of the two wave packets 
can be achieved by shaping the pump field according to 
the dump dipole matrix elements and vice versa: 

E p (cu) cx F d (lu) A (u) 

E d (lu) cx F p (to) A [uj) exp [i<f> D (w)] , (3) 

where A is an arbitrary spectral amplitude, common 
to both fields. Intuitively, this spectral shaping avoids 
pumping of what cannot be dumped (due to a node in 
the dump dipole matrix elements), and vice versa. An ex- 
ample of pump (dump) Franck-Condon factors is shown 
in Fig. |2k(b). The fast oscillation (alternating sign) of 
the Franck-Condon factors for the dump pulse is canceled 
by the relative delay of half a vibration between the two 
pulses, so the pump should only be spectrally shaped ac- 
cording to the slow envelope of Fd(uj). Such shaping is 
easily achieved with current ultrafast pulse shapers |l2| . 
It may seem surprising that a broadband pulse can dump 
a broadband wave-packet to a single state, but since we 
tailored the wave-packet to the pulse, such that all the 
transition amplitudes from levels within the wave-packet 
interfere constructively only at the target state, this is no 
contradiction. While the shaped pump-dump pulse pairs 
achieve spectral selection of a single vibrational level, co- 
herent accumulation is key to further refine the spectral 
resolution to address rotation and hyperfine levels, as well 
as to accomplish near unity transfer efficiency. 



To check the viability of our scheme, we numerically 
simulated the molecular dynamics driven by a train 
of pulses as discussed. The simulation is based on a 
split-operator code 0| that solves the time-dependent 
Schrodinger equation (within the rotating wave approxi- 
mation) with three wave-packets on three potential sur- 
faces, coupled by two arbitrary pulses. According to Eq. 
iJHJ, the pulses can be shaped in two stages: differen- 
tial shaping that ensures overlap of the wave-packets by 
matching the pulses to the different Franck Condon spec- 
tral responses; and common shaping, which affects the 
overall shape of both wave-packets (e.g. common chirp- 
ing). First, the desired effect of differential shaping was 
verified. Indeed, with weak pulses (< n/10) shaped ac- 
cording to the Spectral Franck-Condon function shown 
in Fig. [2 the overlap of the pumped and dumped wave- 
packets was practically unity (> 0.999). We then ex- 
plored the effect of common shaping on the pump-dump 
process. Within the perturbative discussion relevant to 
Eq. J2J, the common spectral amplitude A(u>) is com- 
pletely arbitrary, and for weak enough pulses, this is veri- 
fied by our simulation. Yet, it is desirable to minimize the 
total number of pulses for both practical reasons (more 
tolerant phase locks) and fundamental ones (the total in- 
teraction time is limited by the coherence time of the 
input state), so how strong can the pulses be and still 
qualify as "weak"? The limiting power level is where 
two-photon (Raman) processes by one pump pulse be- 
come pronounced, and here spectral shaping will have an 
important effect. 

In many cases of molecular dynamics, positively 
chirped (red to blue) excitation pulses can strongly sup- 
press Raman processes that adversely affect the input 
wave-packet during the pulse |l4j , leaving the excitation, 
although strong, essentially one-photon. The reason is 
that within the Franck-Condon window, the excited po- 
tential is usually steeper than the ground potential (e.g. 
excitation from the ground vibrational level, localized at 
the zero slope of the ground potential) , so population ini- 
tially excited by the red part of the pulse cannot be later 
de-excited by the blue part because there are no available 
levels to de-excite to. It is clearly shown in simulation 
that when the pulses are positively chirped to be longer 
than the vibration time of both the input and the target 
states, over 50% of population can be selectively trans- 
ferred between two deeply bound vibrational levels with 
one pump-dump pair. Thus, chirping the pulse can im- 
prove considerably the dumping efficiency for the first 
dump pulse that is necessarily strong because it dumps 
to an empty target state. 

The common chirping can help resolve vibrational 
structure around the target state deep in the molecular 
well, yet it cannot resolve rotational / hyperfine struc- 
ture, and certainly not the dense environment around the 
input state. Here, the combination of coherent control 
and coherent accumulation proves powerful - coherent 
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FIG. 3: Simulation results for the coherent accumulation pro- 
cess, (a) the input state population density, (b) the intensity 
temporal profile of the pulses (since the pulses are strongly 
chirped it also represents the shaped power spectrum), (c) 
and (d) are the input and target wave-packet population den- 
sities after one pulse-pair, and (e) and (f) are the correspond- 
ing results after 40 pulses. The energy density was 4.3a;10 -5 
j/cm 2 per pulse for the pump, and varied from 1.7xW~ 2 to 
5a;10~ 5 j/cm 2 for the dump 



control techniques (pulse shaping, chirping) are used to 
achieve a precise state match between specific initial and 
target states, while coherent accumulation allows high 
spectral selectivity and total transfer efficiency. 

We simulated the full coherent accumulation process 
in various scenarios. Figure |3] shows simulation results 
for the interaction of a train of 40 pulse-pairs with the 
molecule at a repetition time of 10 ns, assuming a 30 ns 
radiative lifetime for the excited states. The pump pulse 
area was fixed to 7r/6.6, exciting about 5.7% of the in- 
put population each time, and the dump area was varied 
to match with the accumulation progress. The pulses 
were ~ 10 nm in bandwidth (100 fs transform limited), 
chirped out to ~ 1.5 ps by dispersion of 500,000 fs 2 . The 
pump pulses depleted >90% of the input population, and 
when the inter-pulse phase was tuned to the Raman con- 
dition, 95% of this population reached the target state. 
The purifying nature of the coherent accumulation pro- 
cess is demonstrated by the obtained final wave-packet - 
practically a single state. 

The coherent accumulation process proves to be quite 
robust against intensity fluctuations. Simulations show 
that scaling the intensity of the pump or the dump pulse 
train (or both) by a factor of two leaves the total transfer 
efficiency constant within a few percent. The exact vari- 
ation of dump pulse area according to the accumulated 



population (Eq.(l)) is also not critical. Even if the dump 
area is kept constant, the transfer efficiency is >50% over 
a range of factor of two in intensity. 

Due to the high density of levels near the input state, 
it is inherently impossible to avoid leakage of population 
to nearby levels through two-photon Raman processes, 
which is exactly why CW techniques, such as STIRAP, 
require caution. Although our comb scheme is no excep- 
tion, the deleterious effects of this leakage are suppressed 
for two reasons. First, for a single pulse-pair the leakage 
is diminished by the use of weak, mainly "one photon" 
pulses. Second, assuming the comb is not matched to 
the energy spacing of nearby levels, after a large number 
of pulses N the leakage process resembles an incoherent 
random walk, thus scaling as y/N, whereas the coherent 
depletion of the input state scales linearly as N, causing 
it to dominate. In our simulation, although the input is 
depleted by 90% after 40 pulses, it remains >98% pure. 

Although we considered wave-packet dynamics in one 
excited potential, the scheme is easily generalized to dy- 
namics in a set of coupled potentials (e.g spin-orbit). The 
pulses should then be designed according to the (more 
complex) shape and phase of the transition dipole ma- 
trix elements to the set of excited potentials. In addi- 
tion, since ultrashort pulses are used, the excited wave- 
packets are broad and deeply bound within the excited 
potential, so their dynamics is fast. Consequently, the 
scheme is immune to small perturbations, such as hyper- 
fine interactions, that affect the inter-atomic potential 
near dissociation, and are usually not well known. 

To conclude, the presented scheme is a unique 
and powerful combination of frequency domain control 
(comb) and time-domain control (molecular dynamics). 
As such, it enables performance of coherent control tasks 
with both high efficiency and unprecedented spectral res- 
olution. Due to the use of weak pulses the process can 
be analyzed within a perturbative model, thus opening 
an analytic path to strong field problems that were so far 
accessible only by numerical optimizations. We believe 
the scheme is very general and will find applications in 
many areas of coherent control. Specifically, it can be 
used to produce deeply bound ultracold molecules. 
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